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We examine cross sections for the processes ^ NKK in the framework of a phenomenological 
Lagrangian. We include contributions from A and E resonances up to spin 3/2, as well as those from 
an exotic O"^. We allow the Q"*" to have spin 1/2 or 3/2, with either positive or negative parity in 
each case. We also allow the state to be either isovector or isoscalar. We find that the scenario that 
' ' most closely matches observations at Jefferson Laboratory requires a moderately large coupling of 
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I. INTRODUCTION 

(N ; 

^ , A. Experimental Indications 
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. In the past several months, a number of experimental groups have reported signals for a pentaquark state called 
00 ' the 0+ [1] - [9]. The first evidence for such a state was reported by the Spring-8 Collaboration [1]. The search by 
, Spring-8 was motivated by predictions, made within the framework of the chiral soliton model, by Diakonov et al. 
'nI • [10]. Most searches that have reported evidence for the state put its mass around 1540 MeV. However, in all cases, 
] the experimental resolution has been such that only upper limits for the width of the state could be given. Evidence 
for other pentaquarks predicted as partners to the Q'^, particularly the S , has also been reported by the NA49 
^ \ Collaboration [11]. Using the time-delay technique, Kelkar et al. have found evidence of not only the G"*" in K^N 

■ scattering data, but also a possible spin-orbit partner [12] along with a third possible state. 
^ ' Despite the number of pentaquark sightings, the situation is far from clear. Two members of the NA49 Collaboration 
(~| , produced a minority report pointing out that there was no strong evidence for the existence of the S in older, 
higher-precision data [13]. The HERA-B Collaboration sees no evidence for the 8+ [14], and the BES Collaboration 
also report no evidence in their searches [15]. Searches at RHIC have also yielded no evidence so far [16]. Because 
some experiments have reported signals for the pentaquarks, while others have seen no sign of them, Karliner and 



X 



5h ' Lipkin [17] have postulated the existence of a 'crypoexotic' N* that plays a significant role in production of the Q~^. 
. - t In addition, none of the experiments that report a signal for any of the pentaquarks can say anything about their 
spin or parity. 

Quite apart from the question of the existence of the 8"*", the question of its width is also very interesting. Nussinov 
[18] has examined the implication of such a state for existing K~^d data, and has concluded that the width of the state 
had to be less than 6 MeV. Arndt and collaborators [19] have performed a similar analysis on K^N scattering data, 
and have concluded that the width has to be less than 1 MeV, while Haidenbauer and Krein [20] conclude that the 
width of the state must be less than 5 MeV, or that its mass must be much lower than reported. Cahn and Trilling 
[21] have suggested that the width is 0.9 ± 0.3 MeV, based on their analysis of data from A'+ collisions on xenon. 
Gibbs [22] has also examined K^d data and has extracted a width of 0.9 ± 0.2 MeV. 



B. Theoretical Implications 



The existence of a pentaquark state wouldn't be too jarring for most QCD practitioners, as multiquark states 
have been anticipated for decades. However, its light mass and apparently narrow width are difficult to explain in 
a 'conventional' scenario, and have stimulated much discussion and many postulates. Dzierba et al. have raised the 
possibility that the 'signal' is really a kinematic reflection [23]. Jaffe and Wilczek [24] have constructed a diquark 



scenario for the 0"*". One consequence of their scenario is that the state should have a spin-orbit partner, for which 
there is Httle or no evidence to date. Capstick and collaborators [25] have suggested that the state is as narrow as it 
is because it has isospin 2. This means that there should be isospin partners, none of which have been seen. 

Jennings and Maltman [26] have examined pentaquark phenomenology in a number of scenarios, and conclude that 
such a state fits into the quark model picture if its parity is positive, but this implies the existence of spin-orbit 
partners. Karliner and Lipkin [27] invoke the mixing of two nearly-dcgcncratc KN resonances to explain the narrow 
width of the . They have also speculated on the phenomenology of pentaquark states containing a charm quark 
[28] . In addition, there are many papers that examine the phenomenology of pcntaquarks using QCD sum rules [29] , 
various quark models [30] and string theory [31]. A number of unique scenarios have also been proposed [32]. There 
have even been suggestions that the states seen are in fact heptaquarks [33] , or NKtt bound states [34] . Some lattice 
simulations suggest that the parity of the state is negative [35,36], while the work of Chiu and Hsieh suggest that it 
is positive [37]. More recent lattice work reports no signal for the pentaquark state [38]. 

C. Cross Section Ramifications 

Among the many unanswered questions regarding the 9+ (and other pentaquark candidates) is that of the cross 
section for its production in a particular reaction. To address this, a number of authors have examined cross sections for 
producing them in a number of reactions [39]- [48]. Most of these have been aimed at determining the spin and parity of 
the state, but they have all provided estimates for the production cross section. Various photoproductioii mechanisms 
and observables are examined in refs. [39] - [48], while the authors of ref. [49] treat the reaction K^p n^K^n with 
kinematics suited to production of the 6+. The authors of ref. [44] also examine production of the 6+ using a pion 
as the incident particle. 

The calculation of ref. [48] is closest in spirit to the work that we present here. Those authors use a phenomenological 
Lagrangian to describe the reaction jn — > nK^K~ . In addition to the contribution of the 0+, they also included 
two S~ hyperons: the E~(1197) and the E~(1660). These hyperons provided the backgroimd contribution in their 
calculation. Close and Zhao [40] have emphasized the importance of comparing the cross section for producing the 
pentaquark states with those for producing non-exotic hyperons. 

In this manuscript, we examine the cross section for the process 7iV NKK. We include many contributions, and 
examine all of the channels that are allowed. By including a number of contributions, we are able to understand the 
roles played by non-exotic hyperons, and by the (j){lQ2Q). More specifically, since some of the properties of the non- 
exotic hyperons are known, these can be used to get a handle on how big are the cross sections for their production, 
and for the production of the 0+. 

The results of this calculation are relevant to past, presnt and future searches using photon beams. For the published 
results so far, this means the searches at JLab [3,5], the search by the Saphir Collaboration [4] and, of course, the 
search by the Spring-8 Collaboration [1]. We note, however, that in ref. [5], the process studied is 7p — > n'K^K^K~ , 
not any of the ones discused in this manuscript. While the main focus of our discussion will be the JLab searches, we 
will also comment on the other two searches where appropriate. 

The rest of this article is organized as follows. The next two sections focus on establishing the framework for 
the calculation: the general amplitTidc, kinematics and cross section are disciissed in the next section, and the 
phenomenological Lagrangian terms and most of the coupling constants needed for building the model are presented 
in section III. The diagrams representing the contributions that are included in this calculation are also shown in that 
section. We present our results in section IV, and a summary and outlook in section V. 

II. GENERAL AMPLITUDE, KINEMATICS AND CROSS SECTION 

A. Kinematics and Cross Section 

We begin by describing the kinematics of the process, k is the momentum of the photon, pi is that of the target 
nucleon, p2 is that of the scattered nucleon, and qi and 92 are the kaon momenta. Momentum conservation gives 

k+pi=P2 + q\+q2- (1) 

This means that when we construct the amplitude for the process using all the four-vectors at our disposal, we can 
eliminate one of these from consideration. 
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The total ccnter-of-mass (com) energy of the process is \/s, where s = [k + pif' . We may define a variable t as the 
square of the momentum transferee! to the nucleon, namely, t = (p2 — Pi)'^ , and this is related to the scattering angle 
of the nucleon in the c-o-m frame. 

The differential cross section for this process is described in terms of five kinematic variables. These may be, for 
instance, two Lorentz invariants and three angles. One obvious choice for one of the invariants is s. The choice of 
the other four quantities can be fairly arbitrary, and depends on what information is being presented. One; choice is 
the scattering angle of the nucleon, 9, or equivalently, t. For the other three variables, we can choose for example, 
Sj^j^ = {qi + 52)^ and dU*^— = d(d*^—d^*^—. Here, O^— and are determined in the rest frame of the KK pair, 
relative to a z' axis defined by the direction of motion of the pair of kaons. Another equally valid choice would be 
snk = {P2 + QiY and dfi^^, where the solid angle is defined in the rest frame of the nucleon-kaon pair. 

The differential cross section is 

^,-r9 ^"^,.9 = tA^^-t / |M|^(icos6'd$* ^, (2) 

where 9 is the scattering angle of the KK pair relative to the momentum of the incident photon, in the rest frame of 
the initial photon and nucleon. 



B. General Amplitude 

Our starting point is the construction of the most general form for the transition amplitude for this process. While 
the requirements of Lorentz covariance and gauge invariance delimit the form of the amplitude, we find that there is 
nevertheless quite a bit of freedom in the form chosen. The most general form is 



iM = U{p2)e^O^'U{pi) (3) 

where 

O'* = aip'l; + a2P^ + a^ql + 047** + ^ {a^pt + aeP^ + a-jq^ + 037'*) 

+ h [agp'^ + awp^ + aii< + 0127") + ^1^ {aispt + "'up'^ + a^q'^ + ai67'*) • (4) 

Note that we have no terms in '^i nor '^2, as the initial and final nucleons each satisfy 

i)U{p) =mU{p). (5) 

The amplitude coefficients Cj are all functions of the kinematic variables s, s^;^, 9, O* and or whatever 
combination of kinematic variables is chosen. Their exact dependence on each of these variables will be determined 

by the specific model constructed. 

Gauge invariance of the amplitude requires that k^O^ = 0, which leads to the four relations 

axk ■ pi + a2k ■ p2 + a^qi ■ k — 0, (6) 

04 + ask ■ pi + aek ■ p2 + argi ■ k = 0, (7) 

agk-pi + aiok-p2 + anqi-k = 0, (8) 

ai2 + ctisfc • Pi + auk ■ p2 + a^qi ■ k = 0. (9) 

Note that there is no condition on either ag or aie- 

From these equations, we can eliminate four of the amplitude coefficients, leaving us with twelve independent ones, 
or Lorentz-Dirac structures, to describe the amplitude. One choice would be to eliminate ai, a^, ag, 012, which gives 
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(a2 + a\oi\)p2,xP\v + (as + ciii^Ogi/iiPii. 



+ (as + ai3^i)pii/7/:* 



+ (ae + ai4^i)p2/n7i/ + (a? + ai5^i)9i;i7i' - ^(as + ai6^i)7M7i/| F^" , 
where Pi^" = e'^k" — e^k^^. Another choice is ai, 05, ag, aia, giving 
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(a2 + ^ae + iiaw + i\^aii)p2^pi,, + (04 + ^1012)^11/7^ 



Pi ■ k 

+ (as + l/tar + ^lan + ^1^015) g'l^pi,. 



^(a8 + ai6^i)7/.7>^K^"- 



Note that these two forms contain a potential kinematic singularity at pi ■ k=0. However, this singularity is outside the 
physically accessible region for the process we are discussing, and since this calculation involves no loop integrations, 
such singularities are of no real concern. 

A further four structures may be eliminated by use of so-called equivalence relations, leaving a total of eight. A 
more detailed discussion of this is beyond the scope of this manuscript. 

III. PHENOMENOLOGICAL LAGRANGIANS 

The framework in which wc treat the process ^fN NKK is the phcnomenological Lagrangian. In this approach, 
all particles are treated as point-like. Their structure is accounted for by inclusion of phcnomenological form factors, 
which we discuss in a later subsection. 

A. Ground State Baryons 

We begin with the Lagrangians needed for the electromagnetic vertices of pseudoscalar mesons and ground state 
baryons. We treat nucleons as an isospin doublet, with N = Kaons are also treated as isospin doublets 

f K+\ 

{K = ( ^0 I ) • 7'" and S are treated as isotriplets. 

In what should be a transparent notation, the electromagnetic part of the Lagrangian is (omitting the ©"•" for the 
time being) 

£1 = F (^-^(1 + nh^A^ + + r3C)7^7.f'''^) N 

+ S (-^{1 + nh^A^ + ^{kf + Tskfh^^^F^'^^ S 

- I [K^ (1 + ra) (d^K) ~ (d^K^) (1 + T3) K] A^^ + H.c. (10) 

where /i^A is the E" ^ A transition magnetic moment, /iA is the magnetic moment of the A, and describe 
the anomalous magnetic moments of the nucleon doublet, and the kf^ are the corresponding quantities for the S 
isotriplet. T3 is the isospin operator for the isotriplet. 

The coupling of pseudoscalar mesons to ground state baryons is described by the Lagrangian 

rj is an isosinglet field representing the rj meson. The last three terms of this Lagrangian are obtained by minimal 
substitution in the first three terms. 

B. Vector Mesons 

The vector mesons that enter into our model are K* and 0. The K* is treated as a vector isodoublet field K^, 
completely analogously to the AT, while the (p is represented by a vector isosinglet field (j)^. The Lagrangian in this 
sector is 
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+ [K^ {d^'K) - (d^'K^) K] (j>^ + [K^ (a^TT • r) - (a'^ift))^ . r] K*. (12) 



C. Baryon Resonances 

There are a number of resonances that need to be taken into account in a calculation such as this. Since the 
experimental target is a nucleon, any of the nucleon or A resonances are expected to play a role. For the energy 
range that we consider, and more particularly, for the scope of this calculation, we find that the most salient points 
can be illustrated without any non-strange resonances. Among the hyperons, any number of them can be included, 
but again we limit the scope so that only the lowest few hyperon resonances are taken into account. In either case, 
we do not consider any baryon with spin greater than 3/2. With the scope of the model limited like this, there are 
only a few Lagrangian terms that must be considered in this sector. The non-exotic hyperons that are included are 
listed in table II. 



1. Spin 1/2 

Lagrangian terms needed for spin- 1/2 resonances are 

(-) (-) (-) 

£4 = Ni^NK^^^Y.* ■ Td^K + Ar^^^7^75 (d^'K) A* + N^mK ^g^'K) 6+ 
rriK rriK rriK 

(-) (-) (^) 

+ N^^MK^ . Td^K + Ni^huL^ (^Qi^K) A* + N^mJL^ {^Qt^K) e_ + H.c, (13) 
rriK niK ttik 

where 0± is the field for 6+ with = 1/2^. The first three terms of this Lagrangian correspond to states with 

= 1/2+ , while the last three terms are for = l/2~. In addition, the 0+ part of the Lagrangian written above 
assumes that the state is an isosinglet. For an isotriplet 6+ with — 1/2+, the Lagrangian would become 

C, = iV^^7M750+ • rd^K + 7V^^7^0_ • rd^K. (14) 
rriK rriK 



2. Spin 3/2 

The Lagrangian terms for spin-3/2 resonances are 

(-) (-) (-) 

niK rriK rriK 

(-) (-) (-) 

+ iV^^^75S* • Tdi'K + NiJ^iK^^ fd^K) A* + 77^^75 (d^K) e_„ + H.c, (15) 
rriK rriK mx 

where the /x indices on the A, S and G fields indicate that they are vector-spinor, spin-3/2 fields. In this calculation, 

we use the Rarita-Schwinger version of such fields. The first three terms arc for resonances with positive parity, while 
the last three are for resonances with negative parity. For an isovector Q the Lagrangian terms are 

rriK rriK 



where the 0±^ represents a 9+ with = 3/2^, respectively. 
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TABLE I. Values of Qbb'm obtained using the Goldberger-Treimann relations. 



Coupling 


/m (GeV) 




Qbb'm 


g-rrNNn 


V2 


1.22 


12.8 


Qnsk 


U'itf 


0.34 


3.2 


Qnak 


u:i6 


-0.718 


-6.51 




« 1.2A 


1.22 


10.37 



D. Coupling Constants 

To evaluate the coupling constants of the ground-state baryons to pseudoscalar mesons, we use the extended 
Goldberger-Treimann relations. For the coupling of the baryons B and B' to the pseudoscalar M, the relation is 



gsB'M 



Ga\ 
Gv) 



Mb + Mb' 



B^B' 



2/ 



M 



(17) 



where /m is the meson decay constant for the pseudoscalar meson M. ( §^ ] is obtained from the semileptonic 

K'^v J B^B' 

decay ot B ^ B' or B' — > B. The values of /m, ( %^ ) (taken from the Review of Particle Physics [50]) and 



B^B' 



gBB'M obtained from those relations are shown in table I. 

The decay width of a vector meson into two pseudoscalars is related to the corresponding coupling constant by 



(18) 



where A(a, b, c) is the K /illen function A(a, fe, c) = + + — 2{ab + ac + be). Prom the measured widths and 
branching fractions of the and K* mesons, we find that 



g4,KK = 4.3, gK'K-K = 5.6. 
In a similar way, the width for the process V — > P7 is 



(19) 



(20) 



which leads to 



g<t>ri-y = 4.3, 5,^,r07 = 0.055, gx^'K"^ = 0.35, gK+'-K+^ = 0.22. 



(21) 



For a baryon B with = 1/2+ , the width for the decay into a pseudoscalar meson P and a ground-state baryon 
B' is 



B^B'P 



^^^'^ {Mb + MB' f \{Mb - MB' f - Ml 



167rM|M|, 



{Ml, Ml,, Ml), 



(22) 



while the corresponding width for a baryon with = 1/2 is 

Tb^b'p 



^^^'^ {Mb - Mb-)"" \{Mb + MbO' - M| 



167rM|M|, 



A^ {Ml, Ml,, Ml). 



(23) 



For baryons with = 3/2+, the widths are 



B^B'P = 



B^B'P = 



g%B'P 



1927rM|M2 
9bb'p 



{Mb + Mb')^ - Ml\ Xi {Ml, Ml,, Ml) 
At {Ml,Ml„Ml) 



1927rM|M2 



{Mb + Mb')-MI 



(24) 
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TABLE II. Values of Qynk for non-exotic hyperons appearing in the model. 



y(Mass) 




r (MeV) 


r 


gVNK 


A(1520) 


3 ~ 
2 


16 


0.45 


15.2 


A(1600) 


1 + 
2 


150 


0.2 


1.05 


A(1670) 


1 - 
2 


35 


0.25 


0.32 


A(1690) 


3 - 
2 


60 


0.25 


5.53 


A(1800) 


1 + 
2 


300 


0.35 


0.86 


A(1810) 


1 + 
2 


150 


0.35 


0.71 


A(1890) 


3 + 
2 


100 


0.3 


1.09 


E(1580) 


3 - 
2 


15 


0.45 


1.95 


E(1620) 


1 - 
2 


80 


0.22 


0.52 


E(1660) 


1 + 
2 


100 


0.2 


0.67 


E(1670) 


3 - 
2 


60 


0.1 


3.88 


E(1750) 


1 - 

\t 


90 


0.26 


0.44 


E(1880) 


2 


80 


0.06 


0.19 


E(1940) 


3 - 
2 


220 


0.13 


3.19 


T.VBLE 111. Valu(>s of (/h ya- for diffort-iit spins. pariti(>s and total widtlis of tlic (-). 


J 


P 


r (MeV) 


geNK 


1 

2 


+ 


1 


0.27 


1 


+ 


10 


0.87 


1 

2 




1 


0.16 


i 

2 




10 


0.50 


•i 
2 


+ 


1 


0.61 


■i 




+ 


10 


1.94 


2 




1 


4.35 
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where the first expression is for a positive parity parent baryon. The non-exotic hyperons that are used in this 
calculation, along with their masses, total widths, spins, parities, and their NK branching fractions and coupling 
constants, obtained from eqns. (22) - (24), are shown in table II. 

As mentioned above, we allow the © pentaquarks to have four different combinations of spin and parity in this 
calculation. In addition, since the width of this particle has not yet been ascertained, we also allow different widths. 
Table III shows the values of the coupling constants we obtain for the different spins, parities and total widths of the 
0, assuming that the NK final state saturates its decays. 

Finally, we note that there are a number of coupling constants for which little information is available. Perhaps 
the most important of these in terms of contributions to the cross sections are the couplings of the vector mesons K* 
and (j), particularly those of the (p to the ground-state baryons. The couplings of the two hyperon resonances that 
lie below the NK threshold, namely the E(1385) and the A(1405) are also not known with much certainty. When 
we display and discuss our results, we will comment further on the effects that these coupling constants have on the 
graphs that we show. 

E. Diagrams 



The diagrams that we include in this calculation are shown in figures 1 to 3. In these diagrams, solid lines represent 

baryons. If a solid line is unlabeled, it represents a nuclcon. Dashed lines represent pscudoscalar mesons, with 
unlabeled dashed lines representing kaons. Wavy lines are photons, and dotted lines are vector mesons. Each diagram 
shown actually represents a set of diagrams, as all allowed permutations of external meson and photon legs are taken 
into account. 
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(a) 




> 




(c) 



FIG. 1. 'Born' diagrams: continuous, unlabeled lines are nucleons. Unlabeled dashed lines are kaons and wavy lines are 
photons. 

We include a number of hyperon resonances in this calculation. These are listed in table II. For each of the 

resonances, there is a corresponding set of diagrams of the kind shown in figure 3. 

There are a number of contributions that have been omitted from this calculation. For instance, we have omitted 
all but the ground-state nucleon, and all of the A resonances. In fact, with the information that is available on how 
these states couple to final states with hidden strangeness, we have found that their contributions to the cross section 
are small. We have also neglected couplings to higher moments of any of the hyperon resonances (figure 3 (d)), as well 
as any contributions that would arise from electromagnetic transitions between excited hyperons and their ground 
states. In principle, there is no a priori reason to expect such contributions to be small, but little is known of those 
couplings. Including such contributions would add too many unknown parameters to the model. 



Apart from the photon, none of the states that enter this calculation are elementary particles: they all have 
substructure, and this substructure is reflected in the need to include some kind of form factor at each interaction 
vertex. Indeed, without such form factors, cross sections grow with energy, and the unitarity limit is quickly violated. 

Inclusion of any form factors in a calculation like this must be done in a manner that preserves gauge invariance, 
and a detailed discussion of all of the issues that arise, and all of the methods and prescriptions for preserving gauge 
invariance, are beyond the scope of this manuscript. In this calculation, we adopt the prescription of assigning an 
overall form factor to gauge invariant sets of diagrams. This means, for instance, that all of the diagrams of the type 
shown in fig. 1 have the same form factor as a multiplicative factor. For all of the form factors, we choose the form 



In this expression, pi is usually the momentum of the off-shell particle with mass to^. In this calculation, we make 
the simplification of setting all of the pj to be equal to s, the total energy in the cm frame, squared. X is chosen 
to be 1.8 GeV, as has been used by other authors. In addition, since we apply this form factor to sets of diagrams, 
we choose mi to be the mass of the lightest ofi^-shell particle in a particular set. The exception to this occurs in the 
diagrams of fig. 2 (f) and (h), where nii is chosen to be the mass of the vector meson in the diagram. The value of the 
integer n depends on the spin of baryons in the diagram. If there are only spin-1/2 baryons in the set of diagrams, n 
is chosen to be unity, while for spin-3/2 baryons, n is chosen to be two. 

The form that we have chosen for the form factors, as well as the manner in which we apply them, is simply a 
prescription, and is not meant to be rigorous. We note that the form factors chosen have the desired effect, producing 
cross sections that are roughly of the correct order of magnitude. Without these form factors, calculated cross sections 
are much too large. 



There are six possible channels to be explored, namely 7p pK^K , jp pK^K , ^p nK^K , 771 

nK~^K~, -yn nK^K and 7n pK^K~ . It will be impossible to present results for all of these channels without 
making this manuscript overly long. We therefore choose a few examples to illustrate the main features of the model. 
We note, however, that the first result reported from JLab used a deuteron target, while searches using proton targets 
have been and are being carried out. Examining the cross sections for both kinds of targets is therefore relevant. 



F. Form Factors 



[42,44] 




(25) 



IV. RESULTS 
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FIG. 2. Diagrams containing vector mesons. The dotted lines represent the vector mesons. 




solid line is a nucleon. In diagram (d), the photon couples to the charge of the intermediate resonance: in this model, we 
neglect couplings to any higher moments of the resonance. 
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In the following subsections, we present the results of our model calculation. We begin with the results of the 
full model, including the (p and A(1520). We then exclude these two states to more closely simulate the pentaquark 
searches that have been carried out at JLab, and examine the effects of the spin, parity, width and isospin of the 
pentaquark on the cross section. We also examine the role of the K* in increasing the cross section for production of 
the e+. 

We note that the S(1385) plays very little role in the results we present, as its contribution to the cross section is 
small. The same is true of the N* states that we consider, as their couplings to AK and T,K final states are generally 
small, at least for the ones we have examined. The A(1405), on the other hand, can significantly affect the cross 
section near threshold in the NK mass distributions. If the value for the coupling of this state to the NK channel 
is chosen to be sufficiently large, a sharp shoulder at lower masses arises in the mass distribution of the NK. The 
absence of such a shoulder in the experimental data limits the size of this coupling constant. We use a value of 5.3 
for this constant. 

In reference [3], the process studied was jd pnK^K^ . The A(1520) was identified in the mass distribution of the 
pK~ pair, the (j) in the K^K~ pair, and the 0+ in the riK^ pair. We assume that either one of the initial nucleons 
takes an active part in the scattering process, while the other acts as a spectator. This would mean that the two 
processes contributing to the pnK^K~ final state are 7p pK^K~ and 7n nK^K~ . In addition, this means 
that in the mass distributions observed, only the proton component of the target contributes to the production of 
the A(1520), and only the neutron component contributes to the production of the isoscalar 0+. The two processes 
IP pK^K~ and — > nK^K~ are therefore the focus of much of our discussion. However, searches in other 
channels have been and are being carried out, and some discussion is devoted to those channels as well. 

A. Full Model 

For all of the results that we display, we present da/dm^p where m?- = {pi +pj)'^, and pi is the momentum of the 
ith particle in the final state. Thus, we expect to see strong resonant effects from the d> in the KK subsystem, and 
similar effects arising from the A(1520) in the NK subsystem. We also expect to see weaker resonant effects from the 
other hyperons that are included in the calculation. 

While the coupling of the A(1520) to NK can be determined from the NK partial width of the state, there is no 
simple way of determining the <pNN coupling constants, except by a detailed analysis of (p photoproduction cross 
sections. Indeed, the angular distributions would have to be analyzed in order to determine the relative magnitudes 
and signs of the vector and tensor couplings. Such an analysis is well beyond the scope of this work. In the results 
that we show for the full calculation, we choose = 4 and Cf = 5. The actual values are of no import for the main 
topic of the article. We must also point out that we have not included any diffractive production of the cp. 

The results that we show are for w = 2.5 GeV. We have examined some of the cross sections for smaller values of 
w, and will comment on those results later in the article. 

1. Isoscalar, Spin- 1/2 0"*" 

Figure 4 shows the differential cross section, da/dMj^j^, for an isoscalar 6+ with spin 1/2. The curves in (a) and 
(b) are for the process nK~^K~ , while (c) is for jp pK'^K~ . In (a) and (b), the width of the 6+ is allowed to 

be 1 McV (solid cnirves) or 10 McV (dashed curves). In addition, (a) results from a 0+ with positive parity, while (b) 
corresponds to one with negative parity. Since the isoscalar 0+ is not produced off the proton in this channel, neither 
its parity nor its width affects the curve that results in (c). The structures seen in this plot arise from kinematic 
reflections of the A(1520) and the (p. 

If there were free neutron targets, the curves in (a) suggest that the 6+ would be relatively easy to observe above 
background, modulo detector efficiency, resolution and acceptance issues. However, for deuteron targets, the presence 
of the proton would modify this somewhat. 

Figure 5 shows the same differential cross sections, but as functions of different invariant masses. Figures 5 (a) 
and (b) show the differential cross sections as a function of the mass of the nucleon-antikaon pair, while (c) and (d) 
show it as a function of the mass of the KK pair. In addition, (a) and (c) result from a neutron target, while (b) 
and (d) are from a proton. In the case of the proton target, the A(1520) dominates the cross section in (b), while the 
contribution of the (j) can be seen as the structure at larger values of Mj^-^. The roles are reversed in (d): the (j) gives 
the prominent peak, while the A(1520) provides the 'plateau' at larger invariant mass. For the neutron target ((a) 
and (c)), A's do not contribute to this channel, but the effects of the S resonances included in the calculation can be 
seen in (a). In this case, the bulk of the cross section comes from the (j), as can be seen from (c). 
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FIG. 4. The difFerential cross section da / dM'^j^ as a function of Mnk, for a spin-1/2 O^. The curves in (a) and (b) are for 
the process -yn —^ nK^K~ , while the curve in (c) is for 7p — > pK'^K~ . In (a) and (b), the solid curves arise for a pentaquark 
with a width of 1 MeV, while the dashed curves correspond to a width of 10 MeV. The curves in (a) arise from a pentaquark 
with positive parity, while those in (b) are for a pentaquark of negative parity. 
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FIG. 5. The differential cross section da/dM'^— as a function of Mj^-^ ((a) and (b)), and dcr/dM'^— as a function of Mj^-^ 
((c) and (d)). The curves in (a) and (c) are for the process ■yn — > nK'^K~ , while those in (b) and (d) are for 7p — > pK'^K~ . 
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FIG. 6. The differential cross section da/dMffj^ as a function of Mnk- The curves in (a) and (c) are for the process 
7P — > pK^K , while those in (b) and (d) are for -yp — > nK~^K . (a) and (b) are for a positive parity Q"*" while (c) and (d) are 
for a 0"*" with negative parity. 



Figure 6 shows the differential cross section, da/dM^j^, for the processes pK^K ((a) and (c)) and jp — > 

nK~^K^ ((b) and (d)). The graphs in (a) and (b) assume that the 6+ has positive parity, while those in (c) and (d) 
are for a pentaquark with negative parity. From these curves, particularly those in (b) and (d), it should be clear 
that observing a O"'" signal could be somewhat problematic unless the contributions from the A(1520) and the (j) were 
excluded. We discuss this in a later subsection. 



2. Isoscalar, Spin-3/2 6+ 

Figure 7 shows the differential cross section if the 0+ is assumed to be an isoscalar with spin 3/2. As with the spin- 
1/2 discussion, the curves in (a) and (b) are for the process 771 nK~^K~ , while those in (c) are for jp pK^K~. 
In (a) and (b), the width of the 6+ is allowed to be 1 MeV (solid curves) or 10 MeV (dashed curves). In addition, 
(a) results from a O"*" with positive parity, while (b) corresponds to one with negative parity. It is interesting to note 
that the height of the peak of the Q'^ for the 3/2~ case is comparable to that of the A(1520) (seen in figure 5 (b), 
for instance) . This is not surprising, since the states are almost degenerate, and the height of the peak at resonance 
depends only on kinematics, which would be largely the same for the two resonances. 

Figure 8 shows the same differential cross section, but as functions of different invariant masses. The curves in (a) 
and (b) show the differential cross sections as a function of the mass of the nucleon-antikaon pair, while the curves in 
(c) and (d) show it as a function of the mass of the KK pair. In each case, the upper graph results from a neutron 
target, while the lower graph is from a proton. Unlike the case with spin 1/2, the contribution to the cross section 
of the 0+ is now significant, especially for the negative parity state, and gives rise to the strong kinematic reflections 
seen in (a), and to a lesser extent, in (c) (the dotted curves, for example). 

Figure 9 shows the differential cross section for the processes jp pK°K ((a) and (c)) and jp nK+K ((b) 
and (d)). The graphs on the left assume that the O"*" has positive parity, while those on the right have negative parity. 
From the curves in (c) and (d), it should be clear that detecting a signal for a 0+ with = 3/2~ would be relatively 
easy. 



3. Isovector, Spin-1/2 6+ 

If the 6+ were isovector, there would be a Q~^~^ state that could be seen in K~^p final states, as well as a 0° 
that could be present in nK^ final states. Figures 10 (a) and (c) show the effect of such a state in jn — > nK~^K~, 
while figures 10 (b) and (d) show the effect in 7p — > pK~^K~. The curves in (a) and (b) assume that the O"*" has 
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FIG. 8. The difFcrcntial cross section as functions of M^-^ ((a) and (b)), or M^-^ ((c) and (d)). The curves in (a) and (c) 
are for the process 7n — > nK~^K~ , while those in (b) and (d) are for 7p — > pK~^K~ . 
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FIG. 9. The differential cross section da jdM^j^ as a function of Mnk- The upper graphs are for the process jp — » pK'^K , 
while the lower ones are for 7p — > nK^K . 

positive parity, while those in (c) and (d) assume that it has negative parity. In all cases, the spin is assumed to 

be 1/2. Figures 10 (e) and (g) show the differential cross section for 7p — > pK^K , while (f) and (h) correspond to 

7P — !■ nK^K . (c) and (f) arc for a 0+ of positive parity, while (g) and (h) assume that it has negative parity. The 
curves in (b) and (d) suggest that a signal for a 6++ should be comparable to that for a G"*", whatever the parity of 
the state. 

B. Omitting A(1520), 0(1020) 

It is clear from the graphs shown in the preceding discussion that the A(1520) and the 0(1020) dominate the cross 
section for 7 A'' — > NKK for most channels. To enhance the possibility of isolating a 0+ signal, experimentalists 
impose kinematic cuts to eliminate the bulk of the contribution from these two states. In our case, we will simple 
eliminate all diagrams containing their contributions from the calculation. The curves that result are presented in 
the next two subsections. 

We note that we could also have imposed the same kinematic cuts on the model. The (background) distributions 
that result when we do this are somewhat different from those that we show, but the salient points of the discussion 
are unchanged. 

1. Isoscalar, Spin-1/2 0"*" 

In figure 11 we show the differential cross section that results for a spin-1/2 6+, when the contributions of the 
A(1520) and 0(1020) are omitted from the calculation. In the case of a positive-parity 0+, a signal that may be 
easily identifiable results. In this figure, the smooth background is provided by the non-exotic hyperons included in 
the calculation. The graphs in (a) and (b) are for 771 — > K~^K~n, while (c) is for -yp — > K^K~p. The curves in (a) 
are for a 8+ with = 1/2+, while those in (b) arise from a 8+ with = 1/2^. 

Figures 12 (a) and (b) show the differential cross section as a function of the invariant mass of the NK pair. The 
effects of the non-exotic hyperon resonances included in the calculation can be seen in these curves. Figures 12 (c) 
and (d) show the same differential cross sections as functions of the invariant mass of the KK pair. Since there are 
no resonances left in this channel (in this model), relatively smooth distributions with no prominent features result. 
In this figure, (a) and (c) arc for "/n nK^ , while (b) and (d) are for 7p pK^K^ . The shoulder seen near 
threshold in (b) results from the sub-threshold A(1405). If a larger coupling constant were chosen for this state, this 
structure would be enhanced, while choosing a sufficiently smaller value will make this feature disappear. 



14 



> 




1.4 1.5 1.6 1.7 1.8 1.9 1.4 1.5 1.6 1.7 1 

FIG. 10. The differential cross sections da/dM^j^ as a function of Mnk, for a spin-3/2 O"*", for the processes 7n nK^ 
((a) and (c)), 7p ^ pK+ K' ((b) and (d)), 7p ^ ^7^°:^° ((e) and (g)) and 7/> ^ n/s:+:^ ((f) and (h)). The curves in (a), (b), 
(e) and (f) all arise from a positive-parity 0"*", while those in (c), (d), (g) and (h) all correspond to a ©"*" of negative parity. 
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FIG. 11. The differential cross section da/dM%j^ as a function of Mnk, for a spin-1/2 The curves in (a) and (b) axe 
for the process 771 — » nK~^K~ , while (c) is for 7p — » pK~^K~ . In addition, (a) is for a 0+ of positive parity, while (b) is for 
one of negative parity. 
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FIG. 12. The differential cross section da/dM^— as a function of Mj^-^ ((a) and (b)), and da / dM^— as a function of M^-^ 
((c) and (d)). The curves in (a) and (c) are for the process 'yn — > nK^ K~ , while those in (b) and (d) are for 7p pK'^ K~ . 



In figure 13 we sliow tlie cross sections for the processes 7p pK^K ((a) and (c)) and 7p nK~^K ((b) and 
(d)), in the channel that would show the isoscalar resonance. The effects of the state are clearly seen, and suggest 
that for a pentaquark of positive parity, either channel should provide a clear signal, while for one of negative parity, 
the channel with two neutral kaons is better. 



2. Isoscalar, Spin-3/2 6+ 

Figure 14 shows the differential cross section for a spin-3/2 0+ for the processes "fn nK^K~ ((a) and (b)), and 
IP ~^ pK^K^ ((c))- The curves in (a) assume that the 0"'" has positive parity, while those in (b) assume negative 
parity. In the case of the negative parity state, its contribution completely dominates the cross section. As mentioned 
before, the signal generated by such a state should be comparable to that generated by the A(1520). The positive 
parity 6+ also provides a large signal above the 'background', although it is not as dominant as in the case of negative 
parity. 

Figure 15 shows the same cross sections in terms of different invariant masses. In (a) and (c) the large signal from 
the 0+, particularly from the negative parity version, show up as large kinematic reflections. 

Figure 16 shows the cross sections for pK^T^ ((a) and (c)) and nK'^T^ ((b) and (d)), both assuming a proton 
target. In all cases, both for positive ((a) and (b)) and negative ((c) and (d)) parity, the model indicates that clear, 
easy-to-isolate signals should be obtainable. 



3. Integrated Cross Sections 

In ref [3], in npK~^K~, 212 events are in the peak for the A(1520), and there are 43 events in the peak of 

the 0+. In addition, the Saphir Collaboration [4] report an estimated cross section of 200 nb for production of the 

0+ in the 7p — > nK~^K channel. In this calculation, if we perform a numerical integration around the peak of the 
A(1520), we find that the cross section in the peak in the channel jp pK~^K~ is of the order of 300 nb, with some 
small fraction of this arising from non-resonant contributions. Table IV shows the integrated cross section under the 
peak of the 0+ for different channels in the different scenarios that we have explored. In each case, the integration is 
performed from M — 2r to M + 2r, where M is the mass of the 0+ and T is its width. 

In each case there is some contribution arising from 'continuum' events that lie in the 'right' kinematic regime. 
These continuum contributions represent a larger portion of the reported cross section for l/2~ pentaquarks than for 
1/2+ pentaquarks, for instance. For a l/2~ pentaquark with a width of 10 MeV, approximately half of the reported 18 
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FIG. 13. The differential cross section — as a function of Mnk- The curves in (a) and (c) are for the process 

° NK 

7P pK^K , while those in (b) and (d) are for 7p — » nK~^K . (a) and (b) are for a positive parity 0"*" while (c) and 
(d) are for a @^ with negative parity. 




FIG. 14. The differential cross section da/dM^j^ as a function of Mnk, for a spin-3/2 . The curves in (a) and (b) are 



for the process jn 
negative. 



%K'^K , while (c) is for 7p pK^ K . In (a), the 0"*" has positive parity, while in (b), its parity is 
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FIG. 16. The differential cross section da j dM^j^ as a function of Mmk- The upper graphs are for the process 7p pK'^K , 
while the lower ones are for 7p —> nK'^K . Those on the left are for a G"*" with positive parity, while those on the right assume 
negative parity for the pentaquark. 
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TABLE IV. Total cross sections for production of the S"*", in different scenarios, for the channels in which it can be produced. 
The numbers in the table are obtained from the versions of the model in which the 4> and A(1520) are omitted. The second to 
fifth columns in the table correspond to a ©"*" with a width of 1 MeV, while the sixth to ninth columns correspond to a width 
of 10 MeV. 



Process 




a (nb), re+ 


= 1 MeV 






a (nb), re+ 


= 10 MeV 






1/2+ 


1/2- 


3/2+ 


3/2" 


1/2+ 


1/2- 


3/2+ 


3/2" 




2.6 


1.0 


0.9 


7.7 


25.3 


9.6 


8.9 


73.6 


7P — > nK~^K 


4.3 


2.4 


2.3 


10.0 


44.8 


27.9 


26.0 


110.5 


•yn — > nK~^K~ 


5.6 


1.7 


2.2 


24.0 


54.5 


18.0 


22.8 


229.9 


jn — > plCK~ 


5.9 


1.8 


2.0 


24.0 


56.1 


19.0 


21.2 


225.1 



nb (in 771 — > nK~^K~) arises from such continuum contributions. For a pontaquark with the same quantum numbers 
but a width of 1 MeV, approximately two-thirds of the reported 1.7 nb in the same channel are from the continuum. 

Assuming that the number of events seen is directly proportional to the cross section, modulo questions of detector 
acceptances, resolution and efficiencies, the JLab numbers suggest that the cross section for the 0+ should be of the 
order of 60 nb aroimd its peak in the channel 777, nK'^K~ . 

The numbers in the table indicate that the versions of the model discussed so far are inconsistent with the 0+ 
signal measured at Jefferson Lab, for instance. The estimated cross section inferred for production of the 0+ in the 
process jn nK^K^ is about 60 nb, and for a pentaquark with ~ 1/2+ and a width of 10 MeV, the cross 
section calculated in this model is 55 nb. However, such a large width for the state appears to be in contradiction 
to cross sections observed in other processes [18] - [22]: consistency with such observations would dictate that the 
preferred scenario is for a pentaquark with a width of 1 MeV. In this case, the scenario that most closely matches 
JLab observations is that with a 3/2" pentaquark. However, the results of this calculation suggests that such a state 
should not need kinematic cuts for observation. None of the scenarios with the narrower 0+ match the reported 
Saphir cross section of 200 nb. 

We have examined cross sections for these processes at smaller values of s/s. While the overall cross sections change, 
the relative strengths of various contributions remain similar to what they are at ^/s = 2.5 GeV. In particular, the 
ratio of integrated cross section for the 0+ and the A(1520), remains similar to what it is at = 2.5 GeV. Thus, the 
discrepancy between the results of our model and the signal seen at JLab would remain as difficult at lower energies. 



C. The Role of the K* 



The preceding subsections suggest that, apart from the case of a 0+ with = 1/2^, a signal for this pentaquark 
should be readily observable, especially when the A(1520) and (p are omitted from the calculation. However, there is 
still an inconsistency between what we have shown and what has been observed experimentally at JLab. 

In the results presented, the contributions of the K* mesons have been limited to diagrams in which they couple 
only to ground state hyperons and nucleons. At this point, there are no contributions in which the K* couple to 
excited hyperons, nor to the 0+. It is relatively easy to include such contributions, and in so doing, we can increase 
the cross section for production of the ©+. 

The phenomenological Lagrangian for the coupling of the ©+ to the K* may be written 

C = N f Gf ^^K; + ^^^^1-^7^7- [O^k;) \ e + H.c, (26) 

if the 9+ is assumed to have = 1/2+ . This is the only scenario we discuss. The two coupling constants ^^'^ 
and Gf'^®^ are unknown. In table V we show results for different values of the vector coupling constant, with the 
tensor coupling set to zero. We see that relatively modest values of the vector coupling are sufficient to give a cross 
section of about 60 nb for production of the 9+ in the 771 nK~^K~ channel. However, even that modest value 
for the coupling constant (of about 6) is somewhat larger than values postulated by some authors. For instance. 

Close and Zhao [40] have suggested that (cf *^®"^) w 3. In other models, similar values have been used. With 
QK'Ne+ _ Q^j. model predicts a very large cross of 280 nb for the production of the 6+ in the 7p — > pK'^T^ 
channel, and a cross section of 73 nb in — > nK^K . This mean means that, in this scenario, the contribution of the 
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K* dominates the production of the 8+. These numbers indicate that the reported JLab resuhs are not consistent 
with the estimated Saphir cross section of 200 nh. In addition, in pK^K , the predicted cross section for production 
of the 6+ is comparable to that for production of the A(1520), implying that the signal should be easily observable. 

TABLE V. Total cross sections for production of the O""", in different scenarios, for different processes. Tlie numbers in the 
tabic arc obtained from the versions of the model in which the (p and A(1520) are omitted. All numbers assume that the 0"*" 
has = 1/2+ and a width of 1 MeV. The tensor coupling is set to zero, and four different values of the vector coupling are 
used. All channels in which the 0"*" can be produced as a resonant state are shown. 



Process 


a (nb), Gf*~»+ =2 


a (nb), Gf"-'^ =A 


a {nb), Gf^^ =6 


a (nb), Gf®^ = 8 


7P -» pK°lf 


32.7 


125.9 


282.2 


502.8 


"fp nK^ K 


11.5 


34.6 


73.4 


128.1 


•yn — > nK^ 


10.0 


30.2 


66.1 


117.8 


7n ^ pK^K- 


30.7 


118.6 


269.7 


484.0 



V. SUMMARY AND OUTLOOK 

We have examined the process KKN within the framework of a phenomenological Lagrangian. We have 

examined a number of scenarios for pentaquark production, and have found that the largest production cross section 
occurs for a O"*" with = 3/2~. However, in such a scenario, the cross section for its production is comparable to 
that for production of the A(1520), and kinematic cuts should probably not be needed to enhance the signal. 

If the 0+ has = 1/2+, the cross section for its production is significantly less than that for production of the 
A(1520), by almost two orders of magnitude if its width is of the order of 1 MeV. This means that special mechanisms 
are required to account for the number of events seen in the JLab experiment, relative to the A(1520). One possibility 
is that the coupling of the pentaquark to the K* is large, so that the dominant mechanism of production involves the 
K*. However, this then leads to a signal for which kinematic cuts should not be necessary. One can also invoke the 
couplings of a number of A''* resonances, but the conclusion about the size of the signal would remain unchanged. 

The only scenario (that we can think of) that would give the appropriate ratio between the cross section for 
production of the A(1520) and the 0+ is for the production of the A(1520) to be suppressed even further than the 
suppression we have already obtained through the use of form factors. However, this seems unlikely, as the calculated 
cross section for producing this state is of the same order of magnitude as those published by Barber et al. [51]. 

A. Outlook 

This calculation is not without its shortcomings. The most important shortcoming is the fact that a very simple 
prescription has been employed to regulate the high-energy behavior of the model. A more realistic treatment, 
consistent with the requirements of gauge invariance, will have to be implemented before such a calculation is applied 
to other processes in the future. 

There are prospects for measuring a number of final states with two pseudoscalar mesons at JLab and at other 
facilities. In particular, there are on-going analyses of the processes jN — > ttttA'', 7A'' — > KKE, jN — > KttA and 
7A'^ Kir's. The calculation we have presented has been set up in such a way that it may be applied to any of 
these (or other) processes in a relatively straightforward manner. The core of the code was originally generated for 
-yN — > ttttN, and the modifications necessary for -yN KKN were not overly difficult. Thus, we may expect to 
apply the methods used herein to other processes in the not-too-distant future. 
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